Dislocation-Mediated Melting in Superfluid Vortex Lattices 
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We describe thermal melting of the two-dimensional vortex lattice in a rotating superfluid by 
generalizing the Halperin and Nelson theory of dislocation-mediated melting, and derive a melting 
temperature proportional to the renormalized shear modulus of the vortex lattice. The rigid-body 
rotation of the superfluid attenuates the effects of lattice compression on the energy of dislocations 
and hence the melting temperature, while not affecting the shearing. Finally, we discuss dislocations 
and thermal melting in inhomogeneous rapidly rotating Bose-Einstein condensates; we delineate 
a phase diagram in the temperature - rotation rate plane, and infer that the thermal melting 
temperature should lie below the Bose-Einstein transition temperature. 
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I. INTRODUCTION 

Experiments on rapidly rotating Bose-Einstein conden- 
sates (EEC) in harmonic magnetic traps have observed 
triangular lattices of vortices [l|, S) S 111- The lattice 
is expected to melt at sufficiently rapid rotation rates, 
ri, close to the radial trap frequency, uj± [a, H S @1' ^.t 
sufficiently high energy temperature [9|. In the rapidly 
rotating regime, the condensate wave function is a lin- 
ear superposition of single particle states of the lowest 
Landau level (LLL) in the Coriolis force; quantum fluc- 
tuations are predicted [a, H S Q to melt the lattice when 
the number of vortices, A'v, is so large that the number 
of particles per vortex (or ffiling factor), i', falls to ~ 
6-10. While quantum fluctuations increase with decreas- 
ing ly, thermal fluctuations at finite temperature increase 
more rapidly. The ratio of the mean square vortex dis- 
placements due to thermal fluctuations to those due to 
quantum fluctuations is - 10^iy-^^^{T/n)\nN^ [1]. Un- 
der typical experimental conditions, melting of the vortex 
lattice induced by thermal fluctuations is more likely to 
be observed before melting due to quantum fluctuations. 

In this paper, we generalize the Halperin and Nel- 
son (HN) theory of dislocation-mediated melting [10| to 
apply to the thermal melting of vortex lattices. The 
HN theory describes thermal melting of lattices as a 
Kosterlitz-Thouless transition [ll| due to unbinding of 
pairs of elastic dislocations in the lattice. The complica- 
tion in generalizing the theory is that the forces between 
vortices in neutral condensates are long ranged, unlike 
the forces in a crystal lattice. The HN theory has been 
applied to vortex lattices in type-II superconductors [12l | , 
where the forces between vortices are short ran ged . It has 
also been discussed in superfluid helium films ^^ under 
the very limiting assumption that the lattice is incom- 
pressible. 

We consider a dilute gas of atoms trapped in an exter- 
nal potential and rotating about the z-axis, taking a zero 
range s-wave interacation between particles with coupling 
constant, g — AttH as/m, where Og is the s-wave scatter- 
ing length for two bosons; we work in the usual mean 
field approximation for the interaction energy. Further- 
more, we describe the elastic behavior of a vortex lattice 



in a rapidly rotating BEG using the coarse-grained, long 
wavelength elastohydrodynamics description 1J| . In lieu 
of treating the vortices as discrete entities, we express 
the deviation of the vortices from their equilibrium po- 
sitions (the triangular lattice sites) as a continuous dis- 
placement field, e. Although the condensates are three 
dimensional, the dominant motion is in the rotational 
plane, to which we confine our attention. The equations 
of motion couple the local (coarse-grained) velocity field, 
V, the local (coarse-grained) particle density, and the vor- 
tex displacement field. In Refs. [Z,^], we used the elasto- 
hydrodynanric equations of motion to derive the modes 
and various correlation functions of system. Here, we use 
the time-independent form of the same equations to un- 
derstand the elastic behavior and thermal melting of a 
superfiuid vortex lattice. 

Dislocations are expected to form through nonlinear 
instabilities of vortex lattice excitations, i.e., Tkachenko 
modes. However, to determine the melting, we need con- 
sider only the thermodynamics of dislocations. In doing 
so, we separate the elastic displacement field into two 
parts: a fiuctuating part from the Tkachenko modes and 
a singular part from the dislocations. In the absence of 
nonlinearities these two terms are independent. 

Prior to deriving the thermal melting of a superfiuid 
vortex lattice, we first review the elastic equations for a 
superfiuid vortex lattice in Sec. |TT1 Then in Sec. IIIIl we 
find the elastic displacement field and energy of a single 
dislocation, and the interaction energy of a dislocation 
pair. In Sec. IIVI we generalize the HN theory of two di- 
mensional melting to a superfluid vortex lattice. Finally, 
in Sec. |Vl we consider melting in an inhomogeneous sys- 
tem. 



II. ELASTICITY 



In this section we discuss the elastic properties of the 
vortex lattice, working in the frame in corotating with the 
lattice. Gonsidering only motion in the transverse plane, 
we write the energy as a function of the two-dimensional 



density n, v and the vortex displacement field, e as [IE 



E{n,v,e} = 



—mnv +nV{r) — —'mnfl r (1) 
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The first term in Eq. ([T]) is the kinetic energy of the super- 
fluid surrounding the vortex lattice, in the second V{r) 
is the trapping potential, and third term is the centrifu- 
gal potential; the fourth term is the interaction energy 
of the BEC, with g2d the usual two-dimensional coupling 
constant (~ Zg, where Z is the axial thickness of the 
condensate); the final term is the elastic energy density. 
To second order in e. 
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The first two terms in Eq. ^ are the compression and 
shear energies, respectively, of the vortex lattice; Ci and 
C2 are the (two-dimensional) density-dependent com- 
pressional and shear moduli (essentially the usual three 
dimensional constants [ij| integrated over the thickness 
of the system in the axial direction) [l^. In the in- 
compressible limit, e.g., in liquid helium, C2 = —C'l = 
nhfl/8 [13, [Hi, while in the LLL hmit (52D" < Ji^), 
C2 ~ 0.119^20'^^ [a, [l3' ^^^ C*! is expected to be small 
[19]. The final term, in which n^ = mfl/TrTi is the 
equilibrium vortex density, is the coupling of the vor- 
tex displacements to density inhomogeneities, with 7 an 
fi-dependent constant [l5|. In the incompressible limit, 
7 = (jrh /to) ln(£/^c), where i = l/y/im^ is the vortex 
(Wigner-Seitz) cell radius, i'^ = h/mH., and ^c is the vor- 
tex core radius [i5[, [20|; in the LLL limit, 7 = nll^/m 

We first calculate the response of the vortex lattice to 
a stationary stress to linear order, assuming that the sys- 
tem is effectively two dimensional with motions only in 
the plane perpendicular to the rotational axis. To linear 
order, thermal and quantum fluctuations of the vortex 
positions do not influence the elastic response. The ve- 
locity field entering Eq. ([1]) is governed, in a stationary 
configuration, by the displacement field of the vortices. 
Conservation of circulation implies that a shift of the po- 
sitions of the vortices from their equilibrium positions 
alters the local velocity of the fluid according to [Tj] : 



V + 2ftxe = — V$, 

TO, 



(3) 



where $ is the non-singular superfluid phase. The curl 
of this equation, (Vxi;)^ = — 2riV • e, relates the trans- 
verse component of the velocity field to the longitudinal 
component of the displacement field. The steady state 
velocity field created by a distortion of the vortex lattice 



obeys V • nv = 0. For a homogeneous system, 'V -v = 0, 
i.e., the longitudinal component of the velocity vanishes. 
Variation of Eq. ([1]) with respect the particle density 
and displacement field gives two first order equations of 
motion, 






n = ^, 



(4) 



where fi is the chemical potential; and 



^] = -4V(CiV-e)-2V-(C2V)e + 7n5;Vn 
= '^ + C, (5) 



where cr is the elastic stress [2l|, and Cn ~ jn'^'Vn is 
the elastic stress due to density inhomogeneity. (Unlike 
in [l5| we separate out the C„ term explicitly in defining 
a.) In general the elastic constants, dependent on the 
density, appear within the derivatives. We first consider 
a homogeneous system, and return to the problem of an 
inhomogeneous system in Sec. El 

The fluid velocity induced by a displacement of the vor- 
tices results in a Coriolis force which balances the elastic 
and external stresses, ^ (including ^„): 



2mnft XV = cr — C,. 



(6) 



The curl and divergence of Eq. ([6]) yield the transverse 
and longitudinal elastic responses, ex and eL, indepen- 
dently. Since V • nv = 0, 



(Tt — 2(72 V £t — Ct, 
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and 



^ai^-2mn^VT = {(^Cx^2C2)'^'^ -^^^nmYi^ = Cl- (8) 
Fourier transforming Eqs. JT]) and ([8]), we find the 



elastic Green's functions, Gij [23], 



GLL(fc) = 
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GLT(fc) = GTL(fc) = 0; 
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transverse stresses give rise only to transverse displace- 
ments, and longitudinal stresses only to longitudinal dis- 
placements. In terms of the elastic Green's functions, the 
dynamic equations for displacement field are [7|, 
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Unlike in conventional elastic systems, the longitudi- 
nal elastic Green's function contains an additional length 
scale, an elastic "penetration depth" : 



2Ci + C2 



(13) 



Note that C2 > while Ci is generally negative [2J] , and 
thus (5^ can be positive or negative. In the incompressible 
regime [2^, S^ = —P/1&, while in the LLL regime, 5"^ « 
Q.mt/(^ = {)m{g2T,n/Kl)e, where (^ = h^/2g2Dnm 
is the healing length in the condensate. In the LLL 
regime, UQ ^ 920^? and hence in both regimes \d\ <C i. 
Thus within the regime of validity of the coarse-grained 
continuum approximation, GllI^) — — l/4J7^n7Ti, and 
GLL(r, r') « {~l/4n^nm)6{r - r'). 

The elastic displacement due to the (longitudinal) driv- 
ing force ^„ = ^n^'^n is thus 



e(r) = I GLL{r,r'Xn{r')d^r' 



' 4D?m 



Vlnn(r), 



(14) 



consistent with the results of Refs. jlSl . |20| . |24| . Correc- 
tions to this result using the elastic Green's functions for 
an inhomogeneous system are of order l/N^. 



III. DISLOCATIONS 



elastic dislocations, which are topological defects in the 
vortex lattice. Such dislocations obey the condition, 



de = b, 



(15) 



for any closed path around the dislocation, where b is 
the Burger's vector [22, |26|. The minimum length of a 
Burger's vector is the lattice spacing. 

We first solve Eqs. ([3]) and dH) for an elastic disloca- 
tion in a superfluid vortex lattice, leaving mathematical 
details for the Appendix. For b — bx and 5 = 0, the 
displacement field is, 
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(16) 



This displacement field is transverse, with vanishing di- 
vergence away from the origin, as expected since the lon- 
gitudinal Green's function vanishes away from the dislo- 
cation. 



To describe dislocation-mediated melting of a super- 
fluid vortex lattice, we first investigate the structure of 



In addition we find that the energy of a single disloca- 
tion is 
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where Kn are modified Bessel functions of the second 
kind, and Cc is the radius of the dislocation core; Cc > 
a, the lattice spacing. We include the dislocation core 
energy in the chemical potential of the system. In the 
limit R ^ 00, the logarithmic term in Eq. ([TT]) is the 



dominant term. 

To second order in the displacement field the interac- 
tions between dislocations are pairwise. The energy of 
two separated dislocations at positions r^ and r^j/, with 
Burger's vectors ba,bai, is, as derived in the Appendix, 
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where Sij is the two dimensional anti-symmetric tensor, The first term in Eq. (jlSp grows with increasing separa- 



tion of the dislocation pairs and is therefore long ranged; 
the final term decreases with the separation at a rate 
depending on 6. In the limit 6 — > oo, where one recov- 
ers classical elasticity, the final term extend to long dis- 
tances; in this limit, Eq. (1181) reduces to the dislocation 
interaction energy of Ref . [io| . 

Two dislocations with opposite Burger's vectors are 
attracted with an interaction which grows logarithmi- 
cally. Dislocations thus bind in pairs at low temperature. 
Dislocation-mediated melting is driven, as we discuss in 
the following section, by unbinding of such pairs. 

The core radius, Cc, the length scale at which the dis- 
location can no longer be described by a continuous dis- 
placement field, as well as the dislocation core energy, can 
only be computed from the local microscopic condensate 
wave function. In the incompressible limit, the energy of 
the vortex lattice depends entirely on the kinetic energy 
of the surrounding superfluid, and thus the core energy is 
linearly proportional to the density. Similarly, in the ro- 
tating frame in the LLL limit, only the interaction energy 
changes with vortex position, and thus the core energy 
is proportional to the square of the density. We approx- 
imate the dislocation core energy by. 



Er. 



ain + a2n 



(19) 



where ai and 02 are positive functions of the rotation 
rate, with the approximate dependence, oi oc fiD.a^ and 
02 oc ga^. Since Ec depends on the position of the vor- 
tices within the dislocation core, where the continuum 
limit is not valid, exact values of ai and 02 must be cal- 
culated in terms of the local condensate wave function. 



IV. DISLOCATION-MEDIATED MELTING 

The dislocation-mediated meltin g of the vortex lattice 
is a Kosterlitz-Thouless transition lli| arising from un- 
binding of pairs of elastic dislocations. Thermal fluctu- 
ations cause dislocations to enter large systems in pairs 
with opposite Burger's vectors. Below the melting tem- 
perature, pairs are bound; pairs of Burger's vectors of 
equal magnitude but opposite direction produce no long 
range elastic deformation, and are thus most stable than 
isolated dislocations. If two opposite dislocations over- 
lap, they mutually annihilate. 

Above the melting temperature, dislocations are un- 
bound, and to a first approximation can be treated in- 
dependently. In the thermodynamic limit (i? — > 00), the 
energy (|17|1 is minimized for a dislocation with the mini- 
mum magnitude of the Burger's vector, the lattice spac- 
ing, a, and thus -Edis — (C2a^/7r) ln(i?/ac) plus a small 
constant. [For a triangular lattice, a? — (2/V3)/nv.] 
This logarithmic form allows us to apply the simple phys- 
ical argument of Kosterlitz and Thouless to determine 
when dislocations are thermodynamically stable in the 
vortex lattice. Since there are ~ {R/acY possible loca- 
tions for a single dislocation in the vortex lattice, the 



entropy of a single dislocation is ^ 21n(_R/ac). Hence, 
the Helmholtz free energy per dislocation is 

F.(^-2T)ln(^). (20) 

Thus for T > C^c? j'li^ = 7;°, dislocations (of size h = 
a) are thermodynamically stable. Isolated dislocations 
remain in the system only above this temperature. 

In the absence of interactions between dislocation 
pairs, T^ would be the melting temperature. As the tem- 
perature increases to T^-, the mean separation between 
two dislocations within a pair di verg es; in general the 
mean separation of a pair a, a' is [11| . 



((r-„-r,0')« / e-'^^-^'^^V^d^ 



(21) 



For T < Tj^ the mean separation is finite, while for 
T > T,°, the logarithmic term in Eq. P^ leads to a 
divergent separation length within a dislocation pair, in- 
dicating that the pair has become unbound. 

To obtain a more accurate description of melting re- 
quires considering interactions of dislocation pairs, which 
can be taken into account by the renormalization-by- 
decimation technique in Refs. [23, |2^. This technique 
requires SjR to be small. We do not give the details 
here, but note that the net effect is to renormalize the 
shear modulus C2 to a temperature dependent function, 
C2.R(r). The shear modulus has an intrinsic dependence 
on the temperature, since it is proportional to the super- 
fluid density. The melting temperature is then [29| 



Tm — 7;— C'2,R(T'm)• 



(22) 



Previous estimates for C2,r for a lattice of particles 
interacting logarithmically (e.g., vortices) have shown 
Ca.R ~ AC2 where A < 1 [11. 

The melting temperature depends only on the shear 
modulus. It is instructive to compare the origin of this 
dependence in a system of vortices with the correspond- 
ing result in a traditional elastic system (e.g., an isotropic 
crystal) [l^. The first two terms of Eq. (fT5|) have the 
same form for the interaction energy of dislocations in 
an elastic system. The third term of Eq. (flS)) , containing 
an additional length scale, 6, is particular to the vortex 
system; this term is a short range interaction, ^ 5^ /r'^. 
In smectic liquid crystals near the nematic-to-smectic- 
A transition, the interaction energy of two dislocations 
also contains an additional length scale, a penetration 
depth; when this length is small enough to screen the 
dislocations in the liquid crystal beyond their core, the 
dislocation pairs unbind at T = [30|. In such a liq- 
uid crystal, all the energy channels (splay, twist, and 
bend) are screened by the compression or stretching of 
the smectic layers. In the superfluid vortex lattice, how- 
ever, only compressional effects are screened (by the rigid 
body motion) while the shearing is unaffected. In the vor- 
tex system the unbinding of the dislocations is promoted 
by shearing alone. 



For temperatures approaching T^ from below, the dis- 
placement correlation function is independent of the sys- 
tem parameters [3]: 



lim 
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In the solid phase, the vortex lattice exhibits algebraic 
long-range order. However, above Tm, the displacement 
correlations function decay roughly exponentially with 
distance, as one sees from renormalization group calcu- 
lations [l(| . 

Extended two dimensional Bose systerns, e.g., helium 
films, or atoms in optical lattices [3l|, |32, [33| , under go a 
(conventional) Kosterlitz-Thouless transition [ll|, [SJ, [S^] 
to a superfluid state at Tkt = 7r?i ns/2m, where mUs = 
Ps is the two-dimensional superfluid mass density. When 
rotating at low temperature, such a system should con- 
tain a vortex lattice, which melts via a dislocation me- 
diated Kosterlitz-Thouless transition at a temperature 
Tm = cl^C2,r/'2it. Using the incompressible shear modu- 
lus in Tin one finds T^/Tkt = (l/47r\/3)ns(Tin)/ns(TKT) 
[13, where n,{T^) / n,{TKT) ~ 1 + 6(1 - T,,JTkt)'^^; 
b in ''He films is measured to be ~ 5.5 j32]. Thus 
Tm/TKT ~ 0.26; renormalization of the shear modulus 
should reduce this result at most by a factor of 2 [12 1. 

Cozzini et al. find from numerical computation of the 
T = shear modulus vs. hn/g2Dn - from the incom- 
pressible regime {Tdl/g2T)n « 1) to the LLL regime 
{hn/g2i)ns >> 1) - that with the exception of a small in- 
crease just above the incompressible regime, C2 ^ nhft/S 
[3a |. This result suggests that T^/Tylt decreases with in- 
creasing hfl/g2T)ns- Preliminary calculations for the pa- 
rameters of the ^^Rb experiment of Ref. [33| give b ~ 4.5 
|39ll40f|. indicating Tm/TxT ^ 0.23 prior to renormaliza- 
tion effects. 



V. MELTING IN AN INHOMOGENEOUS 

SYSTEM 

We expect thermal melting in a rapidly rotating BEC 
in a finite inhomogeneous geometry also to be driven by 
formation of dislocations arising from thermal fluctua- 
tions. Dislocations have in fact been observed in the lat- 
tices in these systems; see Fig. 5a of Ref T|, However, de- 
scribing the structure and unbinding of dislocation pairs 
is much more difficult than in a homogeneous system. On 
the one hand, the elastic constants, Ci and C2, depend 
on the local density, and dislocation energies depend on 
position, while the elastic theory contains second order 
terms explicitly dependent on e in addition to the usual 
terms dependent on derivatives of e [19]. Furthermore, 
in a finite system dislocations can be created not only in 
pairs in the bulk of the system, but an isolated disloca- 
tion can also be formed at the edge, which separates from 
its image dislocation and moves into the bulk. Such ad- 
ditional features must be taken into account numerically 



to determine the details of melting, a problem we leave 
for the future [ilj . 

Here we simply determine an upper bound on the melt- 
ing temperature. We assume that the maximum density, 
no, of the BEC occurs at the center of the trap, where 
a dislocation has maximum energy. Furthermore, this 
maximum energy is less than the energy of a dislocation 
in a condensate of the same size but with uniform den- 
sity, uq. As a dislocation approaches the edge of the trap, 
its energy goes to zero. Thus the dislocation energy in 
an inhomogeneous system is bounded above by 

^dis(r-) < max{^dis} < (^2(0)0^/'^) ln(i?TF/ac), (24) 

where C2(0) is the (unrenormalized) shear modulus at 
the center of the trap, and i?TF is the Thomas-Fermi 
radius of the rotating cloud. Using an upper bound to 
the dislocation energy should yield an upper bound on 
the melting temperature, and thus Eq. (|20|) indicates that 
the temperature at which dislocations arc thcrmodynam- 
ically stable is bounded above by T^ = C2(0)a^/27r. In 
fact, dislocations should appear near the edge of the BEC 
for temperatures below T^ and creep towards the center 
as the T approaches T^^, while for T > T^, we expect 
dislocations to occur throughout the system. 

Melting generally occurs below the Bose-Einstein con- 
densation temperature, Tc « 0.94:h{NuJz{ujj_ - fi^))^/^, 
where w^ and ujj_ are the azimuthal and radial angu- 
lar frequencies of the trap respectively. We first esti- 
mate T^/Tc in the limit of a nearly incompressible fluid 
(hfl <ti gn) in trap geometry, and then consider the 
LLL hmit. In the former regime, T^ = C2{0)a'^/2TT < 
hnn,{0)a'^/16TT = /9^(0)?iV8V3m2, where ps{0) the (two 
dimensional) superfluid mass density at the center of the 
trap. The shear modulus depends on the superfluid mass 
density a function of temperature, which for a weakly in- 
teracting Bose gas is approximately the product of the 
mass, the mean thickness, and the three dimensional con- 
densate density, nTF(^ = 0,T). At the center of the trap 
the condensate density, coarse-grained over the vortices, 
is approximately homogeneous and does not feel the ef- 
fect of the condensate edge; at T = 0, 



nTF(r = 0,r = 0) = 
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(25) 

In writing Eq. ([25]) , we approximate the number of par- 
ticles in the condensate by the total particle number N; 
d± ~ h/y/muTL is the radial trap length, and ^a ~ 1.1596 
is the Abrikosov parameter [4^. At finite tempera- 
ture, UTFir = 0,T) = nTF(r = 0,0)(1 - {T/T.f'). 
The thickness, Z, of the condensate at the center of 
the trap is approximately dz{2gbAn'Y-p{{))/hijJz)^''^ , where 
dz — Ti/ ^mujz is the azimuthal trap length. 

Thus we find T-^jT^ as the solution of 
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(26) 



For the experimental parameters of Refs. [3, Q in ^^R-b 
(flg — 4.8 ± 0.05 nm, d±^ — 3.75 /xm, dz = 4.69 /xm) with 
10'' particles and 17 = 0.9wj_, we have T^/Tc < 0.49. For 
N = 10^ this ratio increases to T,„/Tc < 0.78. If Tm/Tc is 
only slightly less than unity, then the thermal melting of 
the vortex lattice will be hard to observe experimentally. 
Higher rotation rates help to lower the melting temper- 
ature, but in most experiments, the number particles in 
the system decreases too, making thermal melting harder 
to resolve. Note that Eq. ^^ predicts that T^jT^ de- 
creases with increasing scattering length, allowing for the 
possibility of inducing thermal melting via a Feshbach 
resonance. 

In the LLL limit, the density at the center of the trap 
is lii, 



Along the axis, the condensate particles are in the low- 
est harmonic oscillator state, so the thickness of the con- 
densate at the center of the trap is approximately dz- 
Also at the center of the trap, the condensate is in a har- 
monic oscillator state along the axis and approximately 
homogeneous in the plane of rotation, so the density as 
a function of temperature is proportional to 1 — (T/Tc)^, 
which we take as the temperature dependence of Ci in 
the LLL regime. The melting temperature is given by 
the solution of 
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Figure [T] shows schematically the phase diagram of a 
rapidly rotating BEC. The y-axis is the temperature di- 
vided by the condensation temperature, and the x-axis 
is the "rotational rapidity," y = ta,nh.~^{fllu!j_) [43|. In 
this figure we use the parameters of Refs. @,y| as above. 
The curve separating the vortex lattice and vortex liquid 
phases is given by Eq. ([^^ at low rotation, (P5|) at high 
rotation rates, and an interpolation between the two lim- 
its; the two kinks in the curve are artifacts of using these 
two limits with the interpolation. We assume that the 
melting leads directly to a vortex liquid, neglecting a pos- 
sible hexatic phase intermediate between the vortex lat- 
tice and liquid which exhibits no long-range order of the 
vortex positions but maintains algebraic long-range order 
of the triangular orientation of the lattice [10[ . Comput- 
ing the phase boundary between such a hexatic phase and 
a vortex liquid requires microscopic details of the dislo- 
cations (e.g., core radius and energy) beyond the scope 
of this paper. 

The vortex lattice is predicted to melt via quantum 
fluctuations into a new strongly correlated regime. De- 
tails of this regime are discussed in Refs. [51, 10, 144|- With 



the density Eq. ((27)) and v « n/n^, we estimate the phase 
boundary for the strongly correlated system to be at 
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(29) 



In Fig.[Tl where we take a critical filling fraction, fc ^ 10 
[aH]) the phase boundary is indicated by the vertical dot- 
ted line, at 2/ ~ 5.4 (corresponding to il/uj± « 0.99996), 
The phase diagram remains to be studied in detail, e.g, 
near the region where the vortex liquid phase boundary 
intercepts the strongly correlated regime, and the behav- 
ior at higher temperatures. 

At low rotation rates {n/uj± ~ 0.9), the melting tem- 
perature may be too close to the condensation tempera- 
ture to discern the melting of the vortex lattice experi- 
mentally. Higher rotation rates provide a better oppor- 
tunity to observe melting since T^ decreases with Q/u;±. 
However, at sufficiently high rotation rates (e.g., at y ~ 5 
in Fig. [T]), it is not clear how to experimentally distin- 
guish thermal melting from quantum melting, since we 
do not have an adequate understanding of the melted 



states at finite temperature. 
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FIG. 1: Sketch of the phase diagram of rapidly rotating BEG. 
The temperature is in units of the condensation temperature, 
Tc\ the "rotational rapidity" is defined as j/ = tanh~ (Q/lo±) 
[431 . To the right of the dotted line, the system at lower 
temperatures enters a strongly correlated regime. 



of the energy of a dislocation on position. This depen- 
dence is needed to determine whether it is energetically 
preferable for dislocations to form at the edge, or by un- 
binding in the bulk of the condensate. 
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APPENDIX 

In this Appendix we derive the velocity and displace- 
ment fields around an elastic dislocation in a homoge- 
neous two-dimensional vortex lattice, in particular, for 
superfluid vortex lattices, a two-dimensional edge dislo- 
cation [26]. The challenge is to find a solution of Eqs. ([3]) 
and © that satisfies Eq. p5|) for the Burger's vector. 
We first recast Eq. P^ as 



dw 



'jk 



VI. CONCLUSION 

The elasticity of a superfiuid vortex lattice is deter- 
mined by the energy of the vortex cells and the global ro- 
tational fluid fiow; the kinetic energy of the " rigid-body" 
rotational velocity stabilizes the system against a nega- 
tive compression modulus, Ci . Because of the rotational 
motion, the effect of elastic compression is significantly 
weakened, and shearing is the dominant degree of free- 
dom for the vortex lattice. Correspondingly, we find that 
the melting temperature of the vortex lattice depends 
only on the (renormalized) shear modulus. 

Among the unanswered questions for future examina- 
tion are, first, to delineate the properties of the dislo- 
cation core, which will require extensive computation to 
unveil. Dislocation cores play a significant role in finite 
sized systems; in experiments to date, core radii may be 
only an order of magnitude smaller than the system size. 
Also, one needs to take into account a possible transition 
to a hexatic phase, and then the transition of the hcxatic 
to a vortex liquid. The latter process involves the un- 
binding of the dislocations into disclination pairs, whose 
description requires the energy and size of the dislocation 
core [XO!]. 

Dislocations in non-uniform condensates present fur- 
ther unsolved problems, analogous to those for vortices 
in inhomogcncous systems, for example, the dependence 



dxi 



bkS(r -r'), 



(30) 



where e.y is the anti-symmetric tensor, and Wij — 
dcj/dxi. Note that e is not continuous within the entire 
plane (i.e., d'^e/dxdy ^ d'^e/dydx); otherwise, Eq. (I5D1) 
is not valid. 

To satisfy the conditions, W -v — Q and cr + 2mn^xv = 
0, we define two dual (gauge) fields, the stream function, 
"0, and a modified Airy stress function, x- Using the 
elastic stress tensor, o-y, defined from Eq. ([5]) such that 
Gj = —daij/dxi, we write, 



d'4' 
' dxi 



SikSjl 



dxkdxi 



2mnfl6ijip- 



(31) 

(32) 



Substituting these expressions into the curl of Eq. ([3]) and 
Eq. ([50]) . we have. 



J7 _o _o mnCl 



0, 



(34) 



where r' is the position of the dislocation. Taking r' = 0, 
we find 
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The Airy stress function contains two basic terms. For 
finite 6, the first term in Eq. ([35]) is long ranged while the 
second term dies off. The stream function is proportional 
to this second term. For imaginary S we replace ^^^0(2;) 
by —-1^0(2;) (where Yq is a modified Bessel function of 
the first kind) since x('") ^-nd ip{r) are real. [For 6 -^ 
00, we find the typical results for an edge dislocation in 
two dimensions (see Ref [l^), however, this limit is not 
physically relevant here.] 



A. Energy of a Single Dislocation 



According to Eq. p3|) . the displacement field, e, along 
any closed path around a dislocation does not return to 
its initial value. Instead, the displacement field at the 
close of the loop is equal to its initial value plus the 
Burger's vector. Thus, the system includes a "branch 
cut" from the dislocation core to the edge of the system, 
as in Fig. [21 the displacement field for a dislocation is 
single valued. (This structure is similar to that of the 
phase, 4>, around a vortex, which changes from to 27r; 
jumping by 27r across the branch cut.) 



j e-[-2CiV{V 



C2 V^e + nmJl x v\ (Tr, 



where Sj is the normal component of the surface element 
of the section of the boundary (including the branch cut) . 
In the absence of external forces, Eq. ^ implies that 
the integrand in the second integral in Eq. (j38p is zero. 
The integrals at the boundary of the system and around 
the core are negligible; the branch cut contribution is 
dominant. Since e = b on the lower part of the branch 
cut, and e = on the upper part (Fig. [2]), we find finally 



^ei = f[2Cibj{V-e) 



(39) 



+ C2 



' dxi 



' dx. 



fe.(v 



ds: 



where the integral is along the lower branch. Computing 
this expression using Eqs. ((3T|) and ([32]) . relating the 
derivatives of e to uij [2l| , and (I35p and ((36l) , we obtain 

Eq. inD. 



B. Energy of Many Interacting Dislocations 



© 



e=0 



e=b 



To obtain a more realistic description of the system, 
requires taking into account the interactions of disloca- 
tions. To compute the energy of a two-dimensional gas 
of interacting dislocations we first substitute Eq. ([3T|) and 
Eq. (l32|) into Eq. ([37| by relating the derivatives of e to 
Cij |21|], and find. 



FIG. 2; Branch cut extending from the core of a dislocation. 



The energy of an elastic displacement in the superfiuid 
vortex lattice is, 
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where t>T is the transverse flow velocity. Integrating by 
parts and using |t;tP = 2e • OxtJ, from Eq. ((21), we have 
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We consider an aggregate of dislocations with Burger's 
vectors, ba{ra), such that the sum over all Burger's vec- 
tors is zero and x(^) ^^^d its derivatives are zero at infinity 
- unlike a single dislocation, a pair of dislocations with 
equal and opposite Burger's vectors is not long ranged. 
Integrating Eq. (|4T|l by parts we obtain, 



E, 



1 r / 2Ci-fC2 4 rrmn. 
2J ^[ I6C1C2 ^ 



V> d^r, (41) 



and substituting in Eq. 
ing, we derive Eq. JT] 



4Ci 
and Eq. (^51) . and integrat- 
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